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LONG-TERM  GOALS 

The  Navy  has  a  long-term  need  for  diver  and  camera  visibility  predictions  in  order  to  improve  mine 
and  expeditionary  warfare  applications.  The  goal  of  this  research  is  to  evaluate  existing  and  emerging 
diver/camera  visibility  models  prior  to  implementation  into  fleet  operations.  This  research  includes  the 
recommendation  of  instrumentation  and  models  that  can  be  incorporated  into  Navy  survey  and  battle 
groups  for  prediction  and  on-scene  determination  of  diver/camera  visibility  operations. 

OBJECTIVES 

The  Distance  Visibility  Algorithm  (DiVA)  is  a  model  proposed  by  HOBILabs  Inc.,  to  predict  diver  and 
camera  visibility  for  mesoscale  targets.  Originally  the  diver/camera  visibility  was  to  be  dependent  on 
measurements  made  with  only  one  instrument  (a-peta;  HOBILabs  Inc.,)  and  a  lighting  function.  The 
model  now  is  more  complex  with  beam  attenuation,  c,  required  as  well  as  target  properties  and  the  a- 
peta  measurement.  The  objectives  of  this  research  is  to  evaluate  the  mesoscale  size  constraint  of 
DiVA,  to  detennine  the  model’s  assumptions,  check  for  internal  consistencies,  and  then  compare 
DiVA  to  more  classical  fonnulations  of  visibility  such  as  contrast  transmittance  theory  (CTT).  This 
research  also  evaluates  potential  improvements  that  can  be  made  using  a  complex  model  of  diver 
visibility  over  the  classical  CTT  approaches. 

APPROACH 

The  evaluation  of  visibility  models  require  an  in-depth  analysis  of  the  model  theory,  a  definition  of  the 
parameters  used,  and  validation  of  the  predicted  visibility  using  in-situ  measurements.  The  theoretical 
approach  was  to  compare  the  Modulation  Transfer  Function  fonnulation  of  contrast  used  by  DiVA 
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with  the  contrast  transmittance  theory  (CTT)  used  by  previous  investigators  (  Duntley,  1963).  The 
theoretical  analysis  between  DiVA  and  CTT  includes:  1)  comparison  of  CTT  and  DiVA  contrast 
definitions  and  their  relationship,  2)  comparison  of  CTT  and  DiVA  maximum  horizontal  visibility 
equations,  3)  investigation  of  the  meaning  of  Kl,  the  radiance  attenuation  parameter  used  by  DiVA, 
and  4)  the  sensitivity  of  the  parameters  required  in  DiVA.  The  DiVA  model  uses  an  analytical 
expression  for  the  light  field  which,  in  theory,  should  improve  visibility  estimates,  but  the  extent  of  this 
improvement  relative  to  the  uncertainties  in  target  size,  range,  and  reflectivity  requires  evaluation.  The 
issue  is  what  additional  information,  if  any,  is  provided  by  the  DiVA  model  and  how  one  could  use  it. 

Field  data  was  collected  to  evaluate  the  “mesoscale”  condition  for  visibility  models  and  included  a 
clear  and  turbid  water  environment.  For  DiVA  evaluation,  measurements  of  the  parameters  needed 
were  made  with  the  a-Peta  and  c-Peta  instruments.  To  compare  with  other  models  and  to  check  for 
instrumentation  capabilities,  the  beam  transmission,  absorption,  and  backscattering  measurements  were 
also  made  with  other  commercially  available  instruments.  Visibility  measurements  were  made  with 
Navy  divers  and  with  digital  imagery  taken  of  black,  as  well  as,  black  and  white  targets.  The  digital 
imagery  and  the  measured  optical  properties  were  used  in  a  Graphical  User  Interface  specifically 
designed  to  compare  CTT  and  DiVA  results  and  to  evaluate  KL  versus  the  attenuation  coefficient 
measured  by  the  a-Peta  (Ka_.eta). 

WORK  COMPLETED 


A  theoretical  comparison  between  CTT  and  DiVA  formulations  for  maximum  diver  visibility  is  near 
completion.  This  work  includes  a  comparison  of  the  definitions  used  between  DiVA  and  the  CTT,  and 
a  comparison  of  KL  used  by  DiVA,  and  c.  These  analyses  were  required  to  detennine  whether  the 
MTF  formulation,  at  a  mesoscale  level,  significantly  improves  visibility  estimates.  The  utility  of 
having  the  radiance  distribution  provided  by  DiVA  was  also  evaluated. 

In-situ  data  was  collected  and  analyzed  from  two  exercises;  a  clear  water  condition  off  of  Pensacola, 
FL,  and  turbid  water  condition  off  of  Gulfport,  MS.  Visibility  measurements  of  black,  as  well  as, 
black  and  white  targets  (plate  and  spheres)  submerged  at  depth  were  made  by  Navy  Reserve  divers 
(from  the  Gauging  Littoral  Optics  Program,  CNMOC  and  ONR).  These  targets  were  selected  to 
reproduce  the  methodology  used  by  HOBILabs  Inc.,  in  DiVA  development.  The  in-situ  data  also 
included  digital  imagery  of  targets  at  increasing  distances  so  that  the  radiance  of  the  target  relative  to 
the  background  could  be  determined.  The  oceanographic  measurements  include  a-Peta,  c-Peta,  ac-9 
(Wetlabs  Inc.,),  CTD,  sky  and  sea-surface  conditions,  and  Secchi  disc  range. 

RESULTS 


Comparison  of  CTT  and  DiVA  contrast  definitions  and  their  relationships  Although  defined 
differently,  there  is  a  straightforward  relationship  between  both  contrast  definitions. 
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where  CoivAiO)  and  Ccr/(0)  are  the  inherent  contrasts,  and  L To,  a, id  and  L Bo  are  the  inherent  radiances 
of  the  target  and  the  background  respectively.  Equating  the  numerators  of  both  contrast  definitions: 

>  )CDiVA  (0)  =  CCTT  (0)^BO 


if  TO  Lbq  , 


it  can  be  shown  that 


2 


and 


CCTT  (O) 


2C„„,,(Q) 

i-c„,,(o) 


CK,,( 0) 


Cgr  (0) 

2  +  Cc7t(o) 


The  above  relationships  for  target  contrast  are  also  valid  for  the  contrast  threshold  (the  threshold 
beyond  which  one  can  no  longer  discern  an  object  against  its  background;  C mcTT  and  C „,DiVA  )'■ 
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It  is  common  in  the  literature  to  see  a  value  of  0.02  for  the  CTT  contrast  threshold.  The  corresponding 
contrast  threshold  for  DiVA  would  therefore  be  equal  to  .02/(2  +  .02)  «  .01 .  These  different  threshold 
values  are  important  for  a  consistent  comparison  between  DiVA  and  CTT  predictions. 


Comparison  of  CTT  and  DiVA  maximum  horizontal  visibility  equations:  The  maximum  horizontal 
visibility,  Rff  ,  given  by  DiVA  (Maffione  2002)  is: 
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where  A  is  the  “lighting  factor”  of  DiVA.  The  numerator  in  this  equation  can  be  shown  to  be  identical 
to  the  numerator  in  the  CTT  theory  for  horizontal  visibility.  Substituting  the  following  relationship 
into  the  numerator 
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the  equation  for  maximum  visibility  becomes: 
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Making  use  of  the  previous  relationships  between  DiVA  and  CTT  contrasts,  the  bracketed  expression 
becomes: 
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After  substitution,  the  DiVA  and  CTT  results  can  then  be  compared  for  horizontal  visibility  as: 

in  max  (  ln  [C CTT  (0)  /  C  mCTT  ]  „  „  J  /  n  max  )  [C CTT  (0 )  /  CmC7T  ] 

\Kvis  /DiVA  =  - A -  aM  l  Kvis  /CTT  =  - 
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It  can  be  seen  that  the  main  difference  between  DiVA  and  CTT  is  the  assumption  that  the  extinction 
coefficient,  c,  needs  to  be  replaced  by  KL,  the  radiance  attenuation  coefficient.  The  above  is  an 
important  result  because  it  shows  that  DiVA  will,  in  general,  predict  greater  maximum  horizontal 
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visibility  than  CTT  since  Kl  <  c.  In  particular,  for  a  black  target,  with  the  target  contrast  equal  to  -1 .0 
and  the  CTT  contrast  threshold  equal  to  -0.02  (i.e.,  a  DiVA  contrast  threshold  of  ~  -0.01): 
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The  above  result  shows  that,  by  defining  target  contrast  according  to  the  Modulation  Transfer  Function 
(MTF)  approach,  DiVA  complicates  the  algebra  and  results  in  intricate  equations  which  can  obscure 
the  underlying  similarities  between  both  diver  visibility  algorithms. 


Our  analysis  shows  that  the  MTF-based  definition  of  contrast  doesn’t  add  new  information  regarding 
contrast  transmittance.  However  the  DiVA  has  been  put  forward  by  HOBILabs  as  a  “mesoscale” 
visibility  algorithm.  Mesoscale  presumably  implies  low  angular  frequencies  relative  to  the  observer  or 
camera.  For  example,  the  angular  frequency  corresponding  to  a  1  meter  wide  circular  target  at  a 
distance  of  3  meters  would  have  a  value  of  3  cycles/radian.  However,  the  same  size  target  at  a  distance 
of  10  meters  would  correspond  to  an  angular  frequency  of  10  cycles/radian.  Historically  the  MTF- 
based  definition  of  contrast  has  been  reserved  to  treat  and  simulate  the  blurring  of  the  higher  spatial 
frequencies  making  up  the  finer  details  within  a  target  (target  detection  versus  identification).  This 
blurring  is  due  to  the  forward  scattering  of  photons  as  they  propagate  from  the  target  to  the  observer  or 
camera.  There  is  nothing  inherently  wrong  with  using  the  MTF-based  definition  of  contrast  for  target 
discernment,  although  its  use  increases  the  complexity  of  the  algebra  without  providing  additional 
insight  into  the  diver  visibility  problem.  However,  it  is  important  for  DiVA  to  detennine  how  KL 
varies  and  its  relationship  to  c  and  Ka.eta  since  these  are  not  equivalent. 

Investigation  of  Kl  :  The  HOBILabs  report  states  that  Kl  is  an  attenuation  coefficient  that  can  be 
equated  directly  or  indirectly  to  that  measured  by  the  a-Beta  (i.e.  the  Ka.eJ.  The  argument  for  the 
DiVA  algorithm  therefore  rests  on  replacing  c  by  KL  in  addition  to  the  MTF-based  definition  of 
contrast.  This  suggests  that  one  critical  aspect  of  DiVA  rests  in  Kl  and  its  measurement  by  the  a-[lcta. 
The  uncertainty  in  determining  Kl  due  to  target  size,  reflectance,  and  shape  are  all  issues.  This 
attenuation  must  include  the  radiance  field  and  the  target/background  changes  that  occur.  The  DiVA 
model  attempted  to  take  this  into  account  by  the  addition  of  the  lighting  factor,  A.  The  utility  of  the 
lighting  factor  is  an  issue  as  demonstrated  above. 


The  DiVA  concept  does  make  some  important  contributions.  Diva  does  try  to  include  the  influence  of 
the  radiance  distribution  on  visibility.  DiVA  defines  a  parameter  Q t  as  Lbo/Et  as  a  function  of  b/c, 
sun  angle,  and  look  angle;  where  ET  is  the  irradiance  incident  on  a  target.  Physically,  ET/n  represents 
the  radiance  emanating  from  a  white  Lambertian  target.  This  radiance  contribution  influence  is  not 
evident  in  other  treatments  since  assumptions  of  asymptotic  radiance  distributions  are  made.  Another 
contribution  is  the  recognition  that  beam  attenuation,  c,  is  not  necessarily  the  proper  coefficient  for 
visibility  but  rather  this  elusive  KL  ,  which  is  a  difficult  parameter  to  measure  as  DiVA  discovered. 


One  good  illustration  for  selection  of  the  appropriate  visibility  attenuation  coefficient  is  for  large 
targets  that  are  close  to  the  observer,  but  in  turbid  water.  The  effective  angular  frequency  of  the  target 
with  respect  to  the  camera  (distance/target  diameter)  is  a  critical  parameter  in  determining  if  c  is 
appropriate.  For  turbid  water  with  targets  of  reasonable  size  (i.e.  low  spatial  frequency),  beam  c  can  be 
shown  to  underestimate  the  visibility  range  for  the  target  (see  discussion  below  for  20”  black  target). 


In-situ  diver  visibility  analysis:  In  order  to  investigate  the  magnitude  of  KL  with  respect  to  c  and  Ka  .eta 
several  interactive  GUIs  were  developed.  Figure  1  is  a  GUI  developed  to  analyze  the  digital  images 
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collected  during  the  diver  operations.  With  the  GUI,  pixel  values  of  the  background,  black  (or  black 
and  white)  20”  target  images  are  extracted,  then  normalized  with  respect  to  the  background  and  plotted 
in  the  center  graph.  This  data  is  used  to  fit  Kl  to  the  data  using  the  functional  fonn  [  1  -e"  ,  ]  where  R 
is  the  range.  The  graph  on  the  right  uses  the  target  and  background  pixel  values  to  derive  the  contrast 
using  the  CTT  and  DiVA  formulation  (blue  and  red  curves  respectively).  The  measured  data  is  lit  to  a 
curve  using  a  “contrast  attenuation  coefficient”  Kc  by  changing  “gamma,  y.”  In  DiVA,  Kl  is  equated 
to  a+  yb,  where  y  describes  the  scattering  contribution  to  the  path  attenuation  coefficient.  Using 
measurements  of  absorption,  attenuation,  and  Ka  .ela  the  y  is  calculated  and  then  compared  with  the 
value  required  to  lit  that  obtained  using  the  images.  Since  there  is  no  additional  information  in  the 
DiVA  formulation,  both  the  CTT  and  DiVA  respond  the  same  but  with  a  different  magnitude  for 
“contrast”  as  described  above. 


Figure  1.  Graphical  Interface  Analysis  Tool  for  analysis  of  digital  imagery 
[Black  target  digital  image  used  to  determine  CTT  and  DiVA  contrasts  for  comparison  of  CTT  and 
DiVA  contrast  formulations;  and  comparison  of  KL ,  c,  and  Ka.eJ 


The  results  using  the  graphical  interface  tool  are  very  informative  with  an  example  from  the  Gulfport 
exercise  shown  above  for  a  black  target.  In  this  example  the  Kl  derived  from  the  radiances  is  equal  to 
that  obtained  from  the  contrast  formulations  (-0.90).  It  is  between  a  and  c.  Note  that  neither  KL  or  Kc 
is  equal  to  the  Ka.eta.  Using  the  relationship  between  Ka.ela  and  a  and  b  gives  a  gamma  of  0. 15,  whereas 
the  gamma  for  the  contrast  is  about  0.45;  the  increase  is  expected  but  the  magnitude  of  the  increase  is 
unknown.  The  value  of  gamma  and  how  it  varies  is  an  area  of  research.  In  other  examples  the  Kc  is 
very  close  to  c.  The  importance  of  the  variability  between  Kl  and  Kc,  c,  and  Ka  .ela  cannot  be 
understated.  This  implies  that  neither  c  or  Ka  .eta  are  always  the  proper  coefficients  for  large  targets. 
Hence  the  variability  of  KL  with  respect  to  c  and  Ka  .eta  (and  therefore  gamma)  is  an  important  aspect  in 
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determining  whether  complex  models  for  visibility  add  significantly  to  existing  diver  visibility 
algorithms.  Currently,  NAVOCEANO  uses  4/c  as  their  standard  algorithm.  For  black  targets,  Vant- 
Davies  Colley  found  that  4.8 /c  fit  their  data  expanding  oceanic  and  riverine  environments.  For  the 
Gulfport,  MS  and  Pensacola,  FL  exercises  values  ranged  from  3.9 'lc  to  6.1/c  (4.4  for  the  example 
above).  For  our  test  areas,  the  visibility  could  have  been  predicted  to  within  about  25%  using  4.8/c 
without  accounting  for  diver  photopic  response,  light  adaptation,  radiance  field,  or  any  other  of  the 
factors  affecting  visibility.  It  has  yet  to  be  demonstrated  that  new  models  significantly  improve  this. 

There  are  places  in  the  DiVA  model  where  some  disparities  arise.  In  running  DiVA  negative  visibility 
does  arise.  These  are  probably  artifacts  of  the  numerical  radiance  model  used  and  the  difficulty  in 
properly  accounting  for  all  of  the  interactions  between  the  target,  the  observer,  and  the  radiance  field. 
For  general  Navy  use  such  results  would  confuse  a  user  about  the  reliability  of  such  a  model. 

IMPACT/APPLICATIONS 

This  work  will  impact  the  selection  of  instrumentation  and  models  that  could  potentially  be 
incorporated  into  survey  ships  and  the  BattleSpace  Profiler  (BSP).  The  use  of  beam  c  versus  Ka  .ela ,  or 
a  combination  of  these,  is  an  important  decision  for  the  Navy.  The  BSP  is  used  by  Mine  Warfare  to 
provide  METOC  (Meteorological  and  Oceanographic)  relevant  parameters  into  MEDAL  (Master 
Environmental  Data  Library)  which  is  the  Tactical  Decision  Aid  used  by  fleet  operations.  The 
selection  of  a  diver  visibility  algorithm  and  proper  parameters  is  critical.  If  a  complex  model  is 
demonstrated  to  work  better,  but  with  too  many  input  parameters  required,  then  the  utility  of  the  model 
is  either  questionable  or  can  be  used  to  direct  other  research  to  obtain  the  needed  parameters. 

TRANSITIONS 

Presently  there  is  no  fonnal  transition;  however,  the  research  results  are  being  used  by  N85  and 
NAVOCEANO  to  detennine  instrumentation  and  models  used  for  diver  visibility. 
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